How sandpiles spill: Sandpile problem in a thick flow regime 
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We obtain an analytical solution of a one-dimensional sandpile problem in a thick flow regime, 
when it can be formulated in terms of linear equations. It is shown that a space periodicity takes place 
during the sandpile evolution even for a sandpile of only one type of particles. Similar periodicity 
was observed previously for many component sandpiles. Space periods are proportional to an input 
flow of particles ro. We find that the surface angle 9 of the pile reaches it's final critical value (Of) 
from lower values only at long times. The deviation (Of — 6) behaves as (t/ro) -1 '' 2 . 

PACS numbers: 83.10.Hh, 83.70.Fn, 83.10Pp, 46.10,+z 



The increasing interest of the physics community in 
the study of granular materials is in one part related 
with the technological applications of such materials, but 
also, because of its unusual properties |Q-@|. Recently, a 
new intriguing phenomena - spontaneous stratification of 
granular mixtures while they are poured in a quasi-two- 
dimensional Hele-Shaw cell - was observed |9]-[l5). After 
the pouring, grains appear to be arranged in periodical 
layers of different species parallel to the pile surface. Sev- 
eral studies were presented in order to explain the strat- 
ification phenomenon; however a simpler question may 
be asked: is it possible to observe something similar to 
this one in the much simpler case of only one type of 
grains? The answer to this question is not a easy task, 
since the sandpile problem is complicated even for such 
a situation. Fortunately, recent studies of granular flows 
made significant progress from the theoretical point of 
view |lq| . More recently, Boutreux, Raphael and P.-G. 
de Gennes jlTj proposed a phenomenological description 
of some special case of granular flows - a so called thick 
flow regime - that admits an analytical treatment. In pa- 
pers 17 1^] this approach was utilized to consider some 
granular flow configurations. 

We shall use the proposed phenomenological equations 
to answer positively the stated above question and to de- 
scribe a total evolution of the sandpile, that will turn to 
be surprisingly non trivial. We will show that the slope 
of the pile approaches its critical value only at long times 
after a rather complicate relaxation process. In fact, we 
shall reconsider a classical sandpile problem using, may 
be, the simplest possible idea without appealing to more 
refined approaches like, for instance, self-organized criti- 
caUtyjnj. 

The authors of paper [JItJ proposed to describe phe- 
nomenologically one-dimensional thick granular flows 
(the flow thickness supposed to be much higher than the 
grain size) by the following equations: 
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Here h(x, t) is the profile of the static part of the ma- 
terial; r(x, t) is the width of a moving granular layer; 
we assume that the flow is from right to left, all rolling 
grains are supposed to move with an equal velocity v. 
The sum of Eqs. (0), d t (r + h) — vd x r = 0, has the 
form of the continuity equation for grains. The right 
hand parts describe the conversion of the static grains 
to rolling ones and vice versa, the meaning of v u is the 
velocity of the uphill fronts, as we shall see later. Usually 
v u > v p3|. Of is a critical angle, or a so called repose 
angle, that is the angle to which the sandpile will evolve. 
The physical reasons to introduce equations (Q) in such 
a form are described in |l^,^8|. These linear equations 
are much simpler for an analytical treatment than the 
previou sly pro posed nonlinear equations for a thin flow 
regime [|lq , po| -p2| . The general solution of equations (|l|) 
may be written immediately in the form |L7[ : 



r(x, t) = u(x + vt) — w(x — v u t) — v u 9ft, 

v + v u 

h(x,t) — w(x — v u t) + v u 0ft, r > 0; (2) 



h(x,t) = const, 



r = 0, 



(1) 



where u(x) and w(x) are arbitrary functions. We shall 
use it to describe phenomenologically the sandpile evolu- 
tion neglecting, as usual |jl7|,[l8), possible near front de- 
viations from the thick flow regime. Let there be a wall 
at x = 0, and grains be pouring permanently from the 
moment t = at this point, so r(x = 0,t = 0) = ro 
is a boundary condition (ro is the thickness of the in- 
come flow). The sandpile is supposed to expand to the 
left, i.e. to x = — oo. There are no particles at first 
moment so r(x < 0, t = 0) = 0, h(x < 0, t = 0) = 
are the initial conditions. Using Eqs. (||) and the 
above initial and boundary conditions one may obtain 
the following solution of Eq. (|l|) with a front mov- 
ing to the left with the velocity v at the time interval 
< t < t\ = (v + v u )r /(vv u 6f) [see also Fig. 1(a)]: 
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FIG. 1. The evolution of the profiles of rolling grains 
r(x) and static grains h(x) at the region — oo < x < 0. 
r(-oo < x < 0,t = 0) = /i(-oo < x < 0,t = 0) = 0. (a) 
< t < ti, the front moves with the velocity v. (6) t = ti 
[see Fig. 2 and Eq. (4)]. (c) ti < t < t2, the front moves with 
some velocity Vd < v < v u , the breaks of the profiles move 
uphill with the velocity v u . (d) t = t<z, the breaks approach 
the wall at x = 0. (e) t% < t < t$, the front proceeds to 
move with the velocity Vd, the break of r(x) moves downhill 
with the velocity v. Note, that the right linear part of r(x) 
is always motion less. After the break of r(x) overtakes the 
front at t = tz, the general configurations (6) - (e) will repeat 
with a lower front velocity. 



r(x,t) = r + 



v + v u 



x eo + vt) 



h(x, t) 



v + v v 



9 f (x + vt) 



(3) 



Q(x) is the Heaviside function (we do not write the mul- 
tiplier <d(x + vt) in the right hand side part of the second 
equation of Eqs. (|^), since, of course, r(x) and h(x) can 
not be negative). Note, that the slope of the static part 
dh/dx = v u 9f/(v + v u ) is less than the critical slope. 
As a result, at the time t\ one has, 



r(x, t) — h(x, t) = r + 



v + v u 



(4) 



for -[(v+v u )/(vv u )]r /9 f < x < (Fig. 1(b)). These 
equations are used as an initial condition to find the so- 
lutions for the time interval t\ < t < t% (the time t3 will 
appear naturally from the solution, and ti is an interme- 
diate time, the meaning of it can be understood from Fig. 
1). We have to add the following boundary conditions: 



r(x = 0,t > 0) = r , 



r x 



h x 



- v d (t - ti),*i <t <t 3 J =0 
-v d {t-t x ),t x <t<t 3 ) =0. (5) 



v + v u r 
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The second and the third equations of Eqs. ([)]) are the 
conditions for the left front of the pile that is supposed 
to move with an (unknown yet) velocity v&. 

From the condition (0) and (a) and Eq. (||) we find the 
value Vd = v u /(v + 2v u ) for t\ < t < t 3 and restore the 
evolution of the pile in this time interval [see Fig.l(b- 
e)]. Note that the profile breaks move uphill with ve- 
locity v u and move downhill with velocity v. Then we 
repeat the described procedure for the next time interval 
i 3 < t < < 5 , etc. 

After these simple but rather tedious calculations we 
obtain a total solution consisting of linear parts, a struc- 
ture of which one can see from Figs. 1 and 2 and following 
equations. Lowest segmented line in Fig. 2 shows the de- 
pendence on time of the front position. Coordinates of 
the segments are: 



x(far, 



<t<t 



2m+l) — 



m(m + 1) (v + v u ) 2 rg 



2 v u \? m 
t, m = 0,l,2,. 



(6) 



(for m = the time interval is < t < t\) with 
V m = mv + (m + l)v u . The particular times shown in 
Fig. 2 are 



m (v + v u )[(m — l)v + (m + l)v u ] r$ 
m(m + 1) (v + v u ) 2 r 



(7) 
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The front coordinates corresponding to times t^m-i at 
which the front velocity changes its value are 



/. \ jv + Vu) r 
x(t 2m -i) = -m — , m 

Vu Of 



1,2. 



(8) 



Note that they are arranged periodically. Two other 
types of lines are shown in Fig. 2: solid lines x = 
v u (t — t2m) and dashed lines x = —v(t — t2 m ), m = 1, 2, . . . 
The lines of the first type depict the uphill movement of 
the breaks of both profiles r{x) and h(x) [Fig. 1(c)]. The 
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lines of the second type show the downhill movement of locity v u and only the break of r(x) moving downhill with 

the break of the profile r(x) with the velocity v - there the velocity v. 
is no any break of h(x) in this situation [see Fig. 1(e)]. 

Thus, in the problem under consideration, there exist Solutions for all regions of Fig. 2, which are connected 

both breaks of r(x) and h(x) moving uphill with the ve- at these lines, look like: 



2 V u V m f V m V «V -m/y 2 m(m + 1 ) ( V + Vu )2 y 

TO = 1,2, 

h(v u (t - t 2m+2 ) < x < v u (t - t 2m ),0) = — -r Q + " r Of [x + — , m = 0,1,2, ... 

2 v u (m + l)(v + v u ) \ V m J 



r (-v(t - t 2m - 2 ),v u (t - t 2m ) < x < 0) = r H — 0/X, 

m(v + v u ) 

m(m + 1) (v + v u ) 2 r vv u \ r 1 v u V m a ( . vv u \ 

—t < x < -v{t - t 2m ),v u (t - t 2m ) = — H — - — ■ -2 fc*/ x + —t , (9) 



(10) 



Here, to = 0. Inequalities in the right hand parts of 
Eqs. © and © define the areas of validity of the solu- 
tions [see Fig. 2]. Eqs. (^|) and ( |lC| ) describe totally the 
sandpile evolution, see Fig. l(b-e). 
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FIG. 2. The areas of different solutions for the pro- 
files of rolling and static grains [see Eqs. (|^) and (jig)]. 
v/t>ii = 0.3. The lower segmented line shows the depen- 
dence of the front coordinate on time. The solid lines 
x — v(t — <2m), m = 1,2,... depict the uphill movement of 
the r(x) and h{x) breaks. The dashed lines x = —v u (t — <2m) 
show the downhill movement of the break of the r(x) profile. 
The points x m are arranged periodically. 



From Eq. ( jl0| ) and Fig. 2, a space periodicity of the 
process is evident: general shapes of the profiles r(x) and 
h(x) are repeated each time when the front moves by the 
distance (v + v u )ro/(v u 0f) to the left. In fact, the pile is 
increased layer by layer, and the expression for the width 
of these layers is the same as the one for the width of 
stratified layers of different fractions in a two-component 
sandpile that was obtained in papers [jlO 13|. One may 
also see from Eq. ( |l0| ) that heights of the pile at times 



t 2m are also periodic in to: 

h(x = 0, t 2m ) = m r , m = 0, 1, 2 . . . (11) 



It follows from Eq. (|T^) that the slope, 8 = dh/dx, 
of the pile will approach its critical value only at infinite 
time: 



S(t 2m -3 <t< t 2m ) =11 — 



m (v + Vu) 



9 f ,to = 1,2. 



(12) 



(for m = 1, < t < t 2 ) [see Fig. 3]. At t 2m -i < t < 
t2rm i m = 1,2, ... there are two different slopes for two 
parts of the profile h(x). For t 2m < t < t 2m+ i,m = 
0,1,2,... all the profile has the same slope. Thus, for 
long times + v u ) 2 / (vVu)]r / 9 f the slope behaves 

as 



2v 9ft 



Of 



(13) 



and 9 relaxes slowly to its final value 9f by a power 
law. At long times, the coordinate of front is obviously 
x = y / 2vr a t/9f, and its velocity tends to — ^vr / (29 ft). 

We presented the solution to the case of zero initial 
conditions - there were no grains at the first moment. 
The case when there is already a pile with the critical 
slope at t = was also considered by us ]23j] . One may 
show that, after some period during which there still ex- 
ists a part of the pile surface with the critical slope, we 
again will come to the considered above situation. Thus, 
at long times one may observe the same behavior for both 
types of initial conditions. 

In summary, we demonstrated in the case of the thick 
flow regime that a space periodicity takes place during a 
sandpile evolution even in the case of the one-component 
pile. The pile spills during a repeating process: grains 
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are piling layer by layer. The thickness of the layers co- 
incides with the thickness of stratified layers at the two- 
component sandpile problem |l^,|l^]. Thus, in the one- 
component pile, we found a precursor of the stratification 
phenomena. 



FIG. 3. The dependence of the relative pile slope on time. 
v/vu — 0.3. At infinity — » Of. The lines are defined for 
< t < <2, t\ < t < ti, tz < t < t§, etc. At times t\ < t < tz, 
<3 < t < ti, etc., the profile of static grain distribution has 
two parts with different slopes. 

The slope of the pile goes to its final critical value after 
a complicated relaxation. At long times 6r — 6<x \J tq /t. 
Formally speaking, we studied only the one-dimensional 
problem, but an ordinary three-dimensional sandpile 
with an axial symmetry (that means that sand is pour- 
ing on an infinite horizontal plane at a single point) may 
be described by the same equations as Eqs. (I) (the 
coordinate x means the distance from the center now). 
Thus, our results stay also valid for such a pile. The 
following questions remain open. How crucially do our 
results depend on the used phenomenological approach? 
Do the sandpile evolution in a thin flow regime differ 
essentially from the described one? 
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